Can sterile neutrinos be the dark matter? 
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We use the Ly-a forest power spectrum measured by the Sloan Digital Sky Survey (SDSS) and 
high-resolution spectroscopy observations in combination with cosmic microwave background and 
galaxy clustering constraints to place limits on a sterile neutrino as a dark matter candidate in 
the warm dark matter (WDM) scenario. Such a neutrino would be created in the early universe 
through mixing with an active neutrino and would suppress structure on scales smaller than its free 
streaming scale. We ran a series of high-resolution hydrodynamic simulations with varying neutrino 
mass to describe the effect of a sterile neutrino on the Ly-a forest power spectrum. We find that 
the mass limit is > 14keV at 95% c.l. (lOkeV at 99.9%), which is nearly an order of magnitude 
tighter constraint than previously published limits and is above the upper limit allowed by X-ray 
constraints, excluding this candidate as dark matter in this model. The corresponding limit for a 
neutrino that decoupled early while in thermal equilibrium is 2.5keV (95 % c.l.). 

PACS numbers: 98.80.Jk, 98.80.Cq 



One of the major unsolved mysteries in cosmology is 
the nature of the dark matter in the universe. Observa- 
tional evidence points towards cold dark matter (CDM), 
for which random velocities are negligible. Two of the 
leading particle physics candidates, the lightest super- 
symmetric partner and axions, both require extensions 
beyond the standard model. At the same time, neutrino 
experiments over the past decade have shown that neu- 
trinos oscillate from one flavor to another, which is only 
possible if they have mass. Current data from atmo- 
spheric and solar neutrino experiments P, 0] are com- 
patible with mixing between the three active neutrino 
families. Perhaps the simplest way to incorporate these 
neutrino phenomena into the standard model is to add 
right-handed neutrinos, just as for other fermions. 

Given this extension of the standard model it is natu- 
ral to ask if these (almost) sterile right-handed neutrinos 
can also explain the dark matter jj]. At least two sterile 
neutrinos are required to explain the origin of neutrino 
mass and existence of different mass mixing scales in so- 
lar and atmospheric neutrinos, so in a model with three 
families of sterile neutrinos a third one can act as dark 
matter 0. Such neutrinos free stream and erase all fluc- 
tuations on scales smaller than the free streaming length. 
This length is roughly proportional to the temperature 
and inversely proportional to the mass of neutrinos. Thus 
if the neutrino mass is sufficiently high, or the tempera- 
ture sufRciently low, then it acts just like CDM and can 
satisfy all of the observational constraints from structure 
formation. Current constraints require the neutrino mass 
to be above l.SkeV HQ. This is below the 5-8keV up- 
per limits from the absence of detection of X-ray photons 
from radiative decays 0) S IS ^3 ■ A massive neutrino in 
the keV range has also been suggested as a possible ex- 
planation for high pulsar velocities 11] and such a model 
can possibly explain baryon asymmetry in the universe 



A sterile neutrino is not completely sterile if it is to 
provide the origin of mass for active neutrinos: it inter- 
acts with active neutrinos and the interaction strength is 
parametrized by the active-sterile mixing angle O, which 
in this model is required to be very small, O < 10~^. In 
this regime sterile neutrinos never reach thermal equilib- 
rium 13 • In general a sterile neutrino decays into active 
ones, but the lifetime can be well above the age of the uni- 
verse over a broad range of masses and mixing angles of 
interest, so it is effectively stable. If the interaction rate is 
energy independent then the momentum distribution of 
sterile neutrinos is simply a reduced version of the distri- 
bution of active neutrinos Q • In practice the interaction 
rate is not constant over the range of masses of interest, 
because at temperatures above the QCD transition more 
interaction channels become available IS ITsL . In this 
paper we use the latest calculation which however 
has only a minor effect relative to the constant interac- 
tion rate, reducing the derived mass limits by about 10% 

For keV masses of interest the corresponding free- 
streaming length is of order a Megaparsec (Mpc) and 
below. Distinguishing between cold and warm dark mat- 
ter thus requires a sensitive probe of linear fluctuations 
on small scales, but nonlinear evolution erases the initial 
conditions on these scales today. Of the current trac- 
ers of density fluctuations the one that is most suitable 
for WDM is the Ly-a forest 0|. It is measured from 
the absorption observed in quasar spectra by neutral hy- 
drogen in the intergalactic medium and has been shown 
to accurately trace the dark matter distribution T^. It 
probes fluctuations down to sub-Mpc scales at rcdshifts 
between 2 and 4, so nonlinear evolution, while not neg- 
ligible, has not erased all of the primordial information. 
Current WDM constraints from the Lya forest 5j do not 
include the latest measurements of the Lya forest from 
the Sloan Digital Sky Survey (SDSS) [13 111. The goal 
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of this letter is to derive new limits by incorporating these 
observational constraints and combining them with a se- 
ries of new hydrodynamic simulations which accurately 
describe the effect of a massive neutrino on the Ly-a for- 
est. 

The linear theory calculations of WDM using CMB- 
FAST [Tsl result in the matter power spectra shown, rel- 
ative to CDM, in the upper left panel of figure ^ We 
plot the ratio of WDM to CDM power for = 6.5, 10, 
14 and 20keV. Instead of the usual 3d power spectrum 
we plot the corresponding Id projection, which is more 
relevant for the comparison to the Id Lya forest obser- 
vations. One can see the suppression of power on scales 
smaller than the free-streaming length, which depends 
on the neutrino mass. While for > lOkeV there is 
hardly any effect for k < 5h/Mpc in 3d (see e.g. relevant 
figures in the corresponding Id power spectrum 

shows more of an effect because small scale modes in 
3d are projected to large scale modes in Id. For exam- 
ple, for nil, — 20keV there is essentially no effect in 3d 
for k < 3h/Mpc and even at fc = 5h/Mpc the power 
suppression is only 2%. SDSS measurements of the flux 
power in Id do not extend above 2h/Mpc and high reso- 
lution spectra are reliable up to 5h/Mpc. So if one were 
interpreting them as measuring 3d power then it would 
be very difficult to detect neutrino masses in this mass 
range. However, the corresponding Id case in figure ^ 
shows a 3% power suppression at 2h/Mpc and 15% at 
5h/Mpc. This rapidly increases with declining mass, so 
that for a 6.5kcV neutrino the power suppression is 15% 
at k=2h/Mpc and a factor of 2 at k=5h/Mpc. 

Nonlinear evolution and hydrodynamic effects further 
modify the linear predictions, which must be addressed 
with simulations. We ran hydrodynamic simulations for 
a series of neutrino masses ranging from 3.4keV to 20keV. 
Many convergence tests and comparisons between differ- 
ent hydrodynamic codes have been performed, which will 
be presented in a separate publication. These tests con- 
firm the accuracy of the original analysis in ,1ES| , which 
was based on a grid of hydro-PM simulations sparsely cal- 
ibrated with hydrodynamic simulations. For the hydro- 
dynamic simulations in this paper we used the Eulerian 
moving frame TVD+PM code described in . The Eu- 
lerian conservation equations are solved in a frame mov- 
ing with the fluid where numerical Mach numbers are 
minimized, allowing thermodynamic variables to be ac- 
curately calculated for both subsonic and supersonic gas. 
Our standard simulations used 20 Mpc/h boxes with 256"^ 
particles for dark matter and 512"^ cells for gas. We used 
10 Mpc/h boxes with equal or twice this resolution to 
test convergence. 

Simulation results are shown in figure ^ for redshifts 
2, 3 and 4 that span the observational range. We have 
adjusted the level of the UV background to match the 
mean absorption as measured from the data. The re- 
sults show that for rrii, = 20keV there are 1-2 per- 
cent effects at k=2h/Mpc at z=4, increasing to 11% for 
nil, = 6.5keV. At k=5h/Mpc the effects are 6% sup- 
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FIG. 1: Ratio of WDM power spectrum relative to CDM 
shown over the relevant observational range. From left to 
right the sterile neutrino masses are 6.5keV, lOkeV, 14keV 
and 20keV. Top left corner shows Id linear power spectrum, 
while the other 3 panels show the ratios from hydrodynamic 
simulations at redshifts 2, 3 and 4. We used concordance 
cosmology with ^Im = 0.28 and Ho = 71km/s/Mpc. Dashed 
vertical line shows the upper limit on k for SDSS data. 



pression for =20keV mass and a factor of 1.5 for 
nil, =6.5keV. These are redshift dependent: while there 
is little differentiating power between models at low red- 
shift, the differences become significantly larger at high 
redshift, where the mean level of absorption is higher and 
the linear power is better preserved. Finally, we note that 
the suppression of small scale power also affects the large 
scale bias of the flux power spectrum, which explains why 
the ratios do not converge to unity on large scales. 

Since the free-streaming scale for these high mass neu- 
trinos is so small one must be careful to have sufficient 
spatial resolution to capture the suppression of power 
on small scales. This is not a trivial requirement: for 
a lOkeV neutrino the scale corresponding to a 50% sup- 
pression of the WDM transfer function relative to CDM 
occurs at 4 times the cell size in our standard simulations 
(1/4 of the Nyquist frequency - for a 20keV neutrino the 
corresponding number is 1/2). Doubling the resolution 
does not change our results significantly (< 20% of the 
size of the WDM effect), so our simulations have con- 
verged over the relevant mass range, but only barely. It 
is important to note that insufficient resolution weakens 
the constraints, since the power suppression is not cap- 
tured on scales below the Nyquist scale corresponding to 
the cell size. This could be an issue for the SPH simu- 
lations used in Q: for example, for their lOkeV sterile 
neutrino mass simulation (corresponding to 2keV ther- 
mal neutrino in their plots) they find essentially no effect 
in the Lya forest on the observable scales, concluding 
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that such masses cannot be probed by the Lya forest, in 
clear contradiction with our results in figure ^ In their 
simulations a 50% suppression of the 10 keV WDM trans- 
fer function relative to CDM happens at the scale of the 
particle separation. This is a factor of 4 worse resolution 
than our cell size for gas, and a factor of 2 worse than 
our particle spacing for dark matter. Thus much of the 
power suppression is missing already in the initial condi- 
tions and the Lagrangian spatial resolution in the SPH 
simulations cannot restore it (additionally, it is unclear 
that the Lagrangian nature of the SPH should help for 
the near-mean or even underdense gas probed by the Lya 
forest at z ~ 4). 

In addition to the Ly-a forest flux power spectrum 
from SDSS ^3 '^^ have added earlier hig h-resolution 
Lya forest constraints in a weak form |23, |23|. When 
testing the robustness of the derived constraints we also 
include the more recent high-resolution Lya forest data 
pi l2i | . While galaxy clustering and CMB data do not 
constrain WDM, they are useful for constraining the re- 
maining cosmological parameters. We use as inputs the 
SDSS galaxy power spectrum psf and CMB power spec- 
trum from WMAP |2y|. Our analysis is based on the 
Monte Carlo Markov Chain (MCMC) method ^ and 
uses CMBFAST [l3| to output both CMB spectra and 
the corresponding matter power spectra P{k). The out- 
put transfer functions are interpolated onto a grid of sim- 
ulations using the matter power spectra rather than the 
neutrino mass, since it is the matter spectrum that is 
most directly related to the observations. 

Our most general cosmological parameter space has 9 
parameters, which are the Hubble constant, matter and 
baryon density, amplitude, slope and running of the pri- 
mordial power spectrum, tensor to scalar ratio, optical 
depth and neutrino mass. Since in most models of in- 
flation tensors and running are expected to be small we 
also explore the constraints when they are set to zero. 

We compare the theoretical Ppik) directly to the mea- 
sured power spectrum. This is particularly important for 
the WDM analysis, where one cannot use the 3d linear 
power spectrum amplitude and slope constraints as given 
in ^3 , since as emphasized there these are valid only in 
the context of standard CDM models without WDM. The 
Lya forest contains several nuisance parameters which we 
are not interested in for the cosmological analysis, so they 
are marginalized over. These include the UV background 
intensity, temperature-density relation of the gas and the 
filtering length (related to Jeans scale [l^). We also in- 
clude a marginalization over several additional physical 
effects, such as fluctuations in the UV background and 
galactic winds [T9ll29l| . 

Applying the standard MCMC analysis to the WDM 
case we find no evidence of WDM: the limit is > 
14keV at 95 % c.l. (lOkeV at 99.9% c.L). The cor- 
responding limit for neutrinos which were in thermal 
equilibrium at a high temperature, when the universe 
had more degrees of freedom, and then decoupled, is 
m„ < 2.5keV at 95% c.l. This constraint is obtained 
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FIG. 2: To the left are the observed SDSS Lya forest flux 
power spectra as a function of redshift from 2.2 (bottom) to 
4.2 (top) in steps of 0.2. To the right are the power spec- 
tra from the high resolution data compiled at redshifts 2.4, 
3.0 and 3.9. For each redshift the thick lines are from the 
best fitted CDM model, while the (generally lower at high k) 
thin lines are for the corresponding WDM model with 6.5keV 
sterile neutrino. The latter is discrepant with both SDSS and 
high resolution data, with most of the distinguishing power 
coming from higher redshifts. 



in our 9 parameter space, but reducing the parameter 
space to the minimal 7 parameters without running and 
tensors does not change the results. The Lya forest data, 
best fitted CDM model and corresponding WDM model 
for nil, — 6.5keV are shown in figure[21 One can see how 
the suppression of power on small scales in WDM makes 
the fit worse. For this figure, where we have not adjusted 
all the other parameters to their best fitted value, the in- 
crease in with WDM is 77 - when the Lya forest model 
parameters and power spectrum amplitude and slope are 
fitted, Ax^ is still 27. Even without high resolution con- 
straints, the poor fit to the SDSS data is apparent, es- 
pecially at higher redshifts. Removing the high resolu- 
tion data only weakens the bounds by 15%. The con- 
verse however is not true: without SDSS the previously 
found constraint (after 10% adjustment for nonthermal 
momentum distribution) is > l.SkeV (95% c.l.) 0. 
This is because within the high resolution data there are 
degeneracies between WDM and many of the nuisance 
parameters such as the temperature of the IGM, UV flux 
and filtering scale. These can be removed by adding the 
large scale flux power spectrum measured by the SDSS 
data. Finally, we note that using the more recent high- 
resolution Lya forest dataj2^ ,2M does not improve the 
limits obtained above. 

Sterile neutrinos that couple to active ones also decay 
and their radiative decays result in photons with energy 
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peaking at close to one half of the neutrino mass, which 
for keV masses can be searched for in X-rays from either 
clusters or from their cumulative contribution in a ran- 
dom direction. Absence of such X-ray emission in the 
Virgo cluster results in an upper limit on the mass of 
8keV 0] while recent reevaluation of X-ray background 
constraints gives an upper limit of 5keV for the value of 
mixing angle that matches the required density of sterile 
neutrinos [T^- These are all below our 99.9% lower Hmit, 
suggesting sterile neutrinos cannot be the dark matter in 
this model. 

Can the bounds presented here be invalidated by some 
additional physical effect in Lya forest that is not in- 
cluded in our model? This is unlikely, but cannot be 
ruled out completely. There are possible physical effects 
that can in principle affect the Lya forest power spectrum 
and while most of them have been shown to be negligi- 
ble or are already part of our standard analysis |29l l3Clj| , 



there remains a possibility that something else will turn 
out to be important. However, it is important to recog- 
nize how successful is the current model in explaining the 
observations. Any potential effects that may be missing 
in the current analysis are constrained by the remarkable 
agreement of the simplest CDM model with the data. It 
seems unlikely that if we lived in a WDM universe its 
signature were erased exactly by some (yet to be discov- 
ered) physical effect. Barring any such cancellations we 
may conclude that the simplest model of sterile neutri- 
nos as the dark matter is ruled out, since the upper limit 
from their decays and the lower limit from their effect on 
large scale structure no longer leaves an open window. 
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